[E(1,1),E(2,1),E(3,1)]=calculate_(0.01);
[E(1,2),E(2,2),E(3,2)]=calculate_(0.005);
p_rough=zeros(3,1);
for i=1:3
   p_rough(i)= log2(abs(E(i,1)))-log2(abs(E(i,2)));
end
fprintf("%s\n",'Error Table')
fprintf("%8s %6s %15s %15s %15s\n",'','h','FD1','FD2','FD3');
fprintf("%8s %6s %15.6e %15.6e %15.6e\n",'','0.01',E(1,1),E(2,1),E(3,1));
fprintf("%8s %6s %15.6e %15.6e %15.6e\n",'','0.005',E(1,2),E(2,2),E(3,2));
fprintf("%s\n",'Rough Accuracy')
fprintf("%8s %15s %15s %15s\n",'','FD1','FD2','FD3');
fprintf("%8s %15.3e %15.3e %15.3e\n",'',p_rough(1),p_rough(2),p_rough(3));
